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We find an analogy between the single skyrmion state in the quantum Hall system and the BCS 
superconducting state and address that the quantum mechanical origin of the skyrmion is electronic 
pairing. The skyrmion phase is found to be unstable for magnetic fields above the critical field 
B C (T) at temperature T, which is well represented by the relation B c (T)/B c (0) « [1 - (T/T c ) 3 ] 1/2 . 



In the quantum Hall systemJH skyrmions, i.e., charged 
spin textures, have been attracting much attention. 
When Barret and his coworkers found that the elec- 
tron spin polarization drops precipitously on either side 
of v = 1 (v = N/Ncf,, where N is the total numbe* 
of electrons and is the Landau level degeneracy 
the ground state was thought to consist of skyrmion 
which was already theoretically studied through the non- 
linear sigma modeH or the Hartree-Fock approach for 
zero temperature.ta Later, the skyrmion was described bst 
a microscopic model without using the sigma model,tj 
however, electron-electron interactions were assumed to 
be unrealistic point interactions. From finite-size cal- 
culations fop. the electron-electron Coulomb interactions, 
Xie and HeH found the lowest energy states, which are 
nearly degenerated but have different total spins, in 
weak Zeeman coupling limit. These low- lying states 
were used to describe skyrmionS|-,which were suggested 
to be symmetry-breaking statesBQ However, it is still 
desirable to know the quantum mechanical origin of 
skyrmions. Here we find that the quamtum Hall sys- 
tem has a pairing effect, similar to superconducting sys- 
tems. Since the BCS theory has been very successful 
in superconductivity^ we may use the same theoretical 
approach to understand the role of electron-electron in- 
teractions in the quantum Hall system. 

In this work, we present a picture that the skyrmion 
state consists of electron-hole pairs, as the BCS super- 
conducting state made of the Cooper pairs. Without 
the nonlinear sigma model, we are able to derive the 
skyrmion wave function near the filling factor v = L. 
which is identical to that of Fertig and his coworkers 
from the eigenstates of the mean-field pairing Hamilto- 
nian. We find that the skyrmion phase is unstable for 
temperatures or magnetic fields above their critical val- 
ues. 

Let us consider a quantum Hall system of ./V electrons 
with the lowest Landau level degeneracy N,/, = N + 1, 
which has the filling factor (y = 1 — e) slightly less than 
one. The system with v = I+e can be also analyzed using 
the well known particle-hole symmetry. The quantum 
Hall sample is assumed to have no impurities, so that the 
skyrmion size is macroscopic in weak Zeeman coupling 
limit. For sufficiently high magnetic fields, the ground 
state of the system will be in a state "J,} with all the 



spins aligned along the same up direction and one hole at 
the origin. When the applied magnetic field is reduced to 
lower the Zeeman energy, ^/J becomes unstable against 
the creation of an electron-hole pair, which is described 
by c m |C m +i,| I^q), where and c mj(J are the creation 
and annihilation operators for the electron in the angular 
momentum I = m and spin a state. Then, the single pair 
states with different m's can be linearly combined to form 
the new ground state or the low-lying excited states of 
the system. 

In analogy with the BCS pairing Hamiltonian, there 
exists an interaction potential V representing pair-pair 
correlations, which are responsible for the rigidity of the 
pair condensate. By noting that electron-electron inter- 
actions are purely the Coulomb interactions in our the- 
ory, the potential energy V in the Hamiltonian can be 
written as 

V = ^ Kn,ri+l,7M« + l C L,l C n+lJ C ™+l,T C «,l' (1) 
ra^n 

where V^ ni;rn2jITl3im4 is the Coulomb interaction matrix 
between the single particle states with I = mi,TO2,TO3, 
and 7714 in the 2-dimensional system. □ As in the case of 
superconductivity, the electron-hole pairings can lead to 
a phase transition, which is called a skyrmionic transi- 
tion (this notion will be clearly discussed later). In this 
transition, the topological charge Q of the state changes 
abruptly from -I to at the critical temperature T c , 
which depends on the strength of the Zeeman coupling. 
Note that although the Coulomb interaction is purely re- 
pulsive, the expectation value of V can be negative, which 
makes the skyrmion phase stable. Using the Hartree- 
Fock decomposition, one can easily show that the pecu- 
liar form of V allows for the scattering between the states 
with different total spins (S z ), but, it conserves L z — S z , 
where L z is the z-component of the total orbital angular 
momentum. The conservation of L z — S z ensures that V 
represents proper interactions for the low-lying states of 
the system, because as found in Ref. ^, there are almost 
degenerate lowest-energy states with the same value of 
L z — S z , which are well separated from higher energy 
states. Similarly, using the particle-hole symmetry, our 
theory can also be extended to the case of v = 1 + e, 
where the wave function has L z + Srt symmetry and Q 
changes from I to at the transition.!] 
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To calculate the physical quantities, we approximate 
all the interactions by their thermal averages except for 
V in the Hamiltonian. The problem is then reduced 
to finding the properties of the reduced Hamiltonian 

H re d = X/m.CT 6 



777 , (J ^777 , u ^777 , (7 



V, where e m ,cr is the self- 
energy for the electron in the state \m,a). In a mean-field 
concept, we ignore the fluctuations about the expecta- 
tion values of cj n _ ) _ 1 ^c rrl j. This suggests that it will be 
useful to express such a product of operators formally as 

( c m+i T Cm -i — ^m) + an d subsequently neglect quan- 
tities which are bilinear in the presumably small fluctua- 
tion term in parentheses. If we follow this procedure with 
our reduced Hamiltonian, we obtain a model Hamiltonian 

H M — ^ y ^m, crC^ (jCjji.cr 

777, a 

- X!( A ™ C m+l,T Cm <i + A ™ C L,l C n+l,1 + A^&tti), (2) 

777 

where b m 's are to be determined self-consistently from 
the relation 67,7 — (c^^ -^c m ^) and A m 's are the order 
parameters defined as 



A777 — ^ ^ ^777,77+1,77,777+1 (^77^1. ) • 



(3) 



As was done by the Bogoliubov transformation in the 
BCS theory, the Hamiltonian in Eq. (|J) can be diagonal- 
ized by a suitable linear transformation which is specified 
by 



t _ t 1 * t 

C 777,l — V m 7777,0 ' U mlm,l 
Ct77+1,T = - u mTm,0 + W m 7 m ,l, 



(4) 



where the coefficients it™ and v m satisfy litmp + l^ml 2 = 1 
and 7 m ,o and 7 m ,i are the new Fermi operators. Using 
these new operators, the model Hamiltonian in Eq. (j^) 
is expressed as 

H M = ~ Y A mKn + Y E m(lm, 7m,0 ~ 7 ro ,l7m,l) 
m. m 

+ Y r ?™(77 t ,7,07m,0 + 7m,l7m,l), (5) 



where i] m = (e m +i,| + e m ,J.)/2, E m = + A^, and 

671 = (£m,i _ e m+i,T)/ 2 - Here E m and £ m satisfy the 
following relations; 
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If we measure the energies with respect to the chemical 
potential lying between the up- and down-spin levels, we 
can ignore the last term in Eq. (j^) because r\ m is negli- 
gibly small, as compared to E m . 

The ground state \^a) of Hm is then written as 



i* G >-n^ii°)=n^ 



VmC 777+l,T 



)|0), (7) 



where |0) is the vacuum state, and the-jground state en- 
ergy which is evaluated analyticallyutD is found to be 
lower than those for the states with A m = 0. Note that 
the ground state is just equivalent to the skyrmion wave- 
function introduced in Ref. [|. This result suggests that 
the quantum mechanical origin of the skyrmion in the 
quantum Hall system is electronic pairing. 

Since the Hamiltonian Hm in Eq. (|^) describes inde- 
pendent excitations and the "f m operators obey Fermi 
statistics, the probability of the excitation in thermal 
equilibrium is given by (ll l:0 j m ,o) = 1 - (7m,i7m,i> = 
f(E m ), where f(E m ) is the usual Fermi function. Using 
this relation and the linear transformation in Eq. (Q), 
Eq. (||) can be rewritten as 

A m = Y Vm - n+1 ' n ' m + lU ^ Vn( -^rn,llm,l ~ 7m,o7m,o) 
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1 ' A 1 

= 2 Y^ V™,n+l,n,m+ljjr tanh ~(3E n , (8) 

7i " 

where (3 — 1/kgT and the prime denotes m ^ n. In the 
first equality in Eq. (||), the expectation values of the off- 
diagonal terms such as 7^ 7 m .i are neglected because 
their contributions are not important in thermal equilib- 
rium. By solving self-consistently the order parameter 
equation, we can obtain A„j as a function of T for each 
magnetic field and find that the order parameters have 
non-zero values only when B < B c and T < T c , where 
B c and T c are the critical magnetic field and critical tem- 
perature, respectively. Since there is an effective positive 
charge in the vicinity of the origin, £ m has in fact the 
dependence of m, where to's are small integer numbers. 
For computational convenience, assuming £ m = £ for all 
values of m, we can determine £ in unit of the typical 
Coulomb energy e 2 /el from the self-consistent equation 
£ = g+E ex tanh(/?£/2), where g = \g* ^bB / '(e 2 / 'el) is the 
Zeeman coupling parameter, I = y/hc/eB± is the mag- 
netic length, and B and B± are the total applied mgnetic 
field and its normal componenet to the 2D plane, respec- 
tively. Here g* and e are the band effective g-factor and 
the dielectric constant of the host material in the quan- 
tum Hall sample. To determine the unknown parameter 
E ex that comes from the exchange interaction, we first 
calculate the critical value of £ from Eq. (8) at T = 
0, which gives zero values for the order parameters A m . 
Using the known critical value (g c = 0.0265) for the Zee- 
man coupling parameterJ3 at which the skyrmion state 
disappears, the value of E ex is determined so that the 
critical value £ c is equivalent to g c + E ex at T = 0. Then, 
this value of E ex is used for the self-consistent solutions 
of Eq. (8) at finite temperatures. In fact, our approach 
of using E ex as a parameter and a constant value for £ m 
change slightly the magnitudes of A777, however, it does 
not affect the main result of our calculations. 

Once A m (T) is calculated, the temperature-dependent 
excitation energy E m and the occupation number f m = 
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f(E m ) can be determined. Then, the specific heat is 
calculated by C = TdS/dT, where S is the electronic 
entropy of a fermion gas, defined as S = — 2ks J2mi(^ ~ 
/ m )ln(l - f m ) + f m \nf m ]. In the 2D electron system 
with the magnetic field perpendicular to the plane, since 
there is no Meissner effect, the system has no latent heat 
and the entropy is continuous at T c . Fig. [j] shows the 
discontinuity of the specific heat at T c , and the size of 
the discontinuity is readily evaluated as follows; 



AC = (C s - C n ) 



T c — 



/(&) [W(£c)] r riA 2 



knT r 



dT 



(9) 



where £ c = £,(T C ), A 2 = J2 m |Am| 2 , and the subscripts s 
and n denote the skyrmion and normal phases, respec- 
tively. Here we point out that for a given nonzero g, the 
size of the discontinuity in the specific heat AC is finite 
for a single skyrmion state, however, AC/C n is negligible 
in the thermodynamic limit N — > oo. Thus, practically, 
the observation of the phase transition will be only pos- 
sible when the Zeeman coupling is very weak so that the 
skyrmion size becomes macroscopic or many skyrmions 
contribute to the discontinuity of the specific heat. In the 
multi-skyrmion case, when we do not take into account 
inter-skyrmion interactions, the discontinuity of the spe- 
cific heat at the filling factor v may be approximated by 
AC[y) re ACN^v — 1|. It should be addressed that our 
skyrmion phase transition will not be the second-order 
phase transition, because no long-range order will exist 
for a finite-size skyrmion. The skyrmion phase transition 
involves the localized spins of electrons and is not related 
to the long-range order, similarly to the Kondo problem, 
where there appears the quasi-bound state between the 
conducticpi-clectron spin and the localized impurity spin 
below T C H 

As already known for zero-temperature systems, there 
is a critical magnetic field for the skyrmion-normal state 
transition.0 Based on our pairing theory in the quan- 
tum Hall system, we can determine the temperature- 
dependent critical magnetic field B C (T), similarly to the 
the Gorter-Casimir phenomcnological theory for super- 
conductors, which gives the relation B C (T) / B C (Q) re 
1 - (T/T c ) 2 . With Bj_ fixed, B C (T) in our system can be 
expressed as 



g c (T) _ B C (T) 
9c(0) B c (0) 



1/2 



(10) 



where /cbT c is calculated to be about 0.24e 2 /eZ for g = 0. 
This formula for B C (T) is found to be in good agree- 
ment with the results obtained directly from the order 
parameter equation in Eq. (||) (see Fig. ||). For GaAs 
samples, T c for g re is approximately 12.0 a/ B± [Tesla] 
K.PJ If fluctuations or impurity effects are considered, 
the critical temperature is expected to be much lowered, 
however, this is beyond the scope of our work. Since 
fluctuations above the mean field is important in two 



dimensional systems, we expect that the sharp specific 
heat jump as shown in Fig. |l| will not take place in real 
systems. The estimation of to what extent B C (T) and 
AC(T) are smeared out by the fluctuations, thus, will 
make our work more explicit as the Ginzburg criterion 
does in the usual Landau theory. Unfortunately, we are 
unable to find any reliable criterion for the validity of our 
mean field theory at this stage, which will be our future 
work. 

In summary, we have presented an analogy between the 
skyrmion state and the superconducting state. The coun- 
terpart of the Cooper pair in superconductivity is found 
to be the pair which consists of a minority-spin elec- 
tron and a majority-spin hole. We find the temperature- 
dependent critical magnetic field B C (T), above which the 
skyrmion is unstable at temperature T. Since our ap- 
proach is based on the mean field approximation and the 
fluctuations above the mean field are usually important 
in two dimensional systems, our calculated B C (T) is con- 
sidered to be an upper bound for the true temperature- 
dependent critical magnetic field. 
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FIG. 1. The relative specific heat C 3 — C n is drawn as a 
function of T for the Zeeman coupling parameter g — 0.01. 

FIG. 2. The temperature-dependent critical magnetic 
field scaled by B c (0) is shown as a function of T/T c , with 
B± fixed. The dotted and solid lines represent the results 
from the order parameter equation and the fitting formula 
1-(T/T C ) 3 , respectively. 
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